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Abstract We consider a randomly perturbed Korteweg-de Vries equation. The perturbation 
is a random potential depending both on space and time, with a white noise behavior in time, 
and a regular, but stationary behavior in space. We investigate the dynamics of the soliton of 
the KdV equation in the presence of this random perturbation, assuming that the amplitude 
of the perturbation is small. We estimate precisely the exit time of the perturbed solution 
from a neighborhood of the modulated soliton, and we obtain the modulation equations for 
the soliton parameters. We moreover prove a central limit theorem for the dispersive part of 
the solution, and investigate the asymptotic behavior in time of the limit process. 

1. Introduction 

Our aim is to describe the dynamics of a soliton solution of the Korteweg-de Vries equation 
in the presence of a random potential, depending both on space and time and which is white 
in time. After the first paper [21] showing "superdiffusion" of the soliton of the KdV equation 
in the presence of an external force which is a white noise in time (see also [1], |16|). the 
interest in such questions of soliton dynamics in the presence of either deterministic or random 
perturbations has recently increased in the mathematical community. In |l5j , e.g. the question 
is investigated with the help of inverse scattering methods, for different types of time-white 
noise perturbations, still for the KdV equation, while in [IT], [12], the case of a soliton of the 
NLS equation is studied, with the presence of a slowly varying deterministic external potential. 
Random potential perturbations for NLS equations have also been considered in [14] and [9]. 
The diffusion of solitons of the KdV equation in the presence of additive noise was numerically 
investigated in [19]. Also, in [5j, we studied the soliton dynamics for a KdV equation with 
an additive space-time noise. Our aim here is to reproduce the analysis of [5] in the case 
of a random potential, which is stationary in space : the solution of the stochastic equation 
starting from a soliton at initial time will then stay close to a modulated soliton up to times 
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small compared to where e is the amplitude of the random pertm'bation (see below). In 
the present case, where the noise is multiplicative (the random potential) we are then able to 
analyze more precisely the modulation equations for the soliton parameters and the linearized 
equation for the remaining (dispersive) part of the solution, and especially its asymptotic 
behavior in time. 

We consider a stochastic KdV equation which may be written in Ito form as 

(1.1) du + {d^u + ^d^{u'^))dt = eudW 

where e > is a small parameter, u is a random process defined on (t, x) £ x M, is 
a Wiener process on L2(M) whose covariance operator 00* is such that is a convolution 
operator on L^(M) defined by 

ct>f{x) = [ k{x- y)f{y)dy, for / G L\R). 
Jr 

The convolution kernel k satisfies 

(1.2) ||fc||i := [ (k^ + ik'f)dx < +00. 

Jr 

Considering a complete orthonormal system (ei)igN in L'^i^), we may alternatively write W 

as 

(1.3) W{t,x) = Y,(3i{t)ct>e^{x), 

(/3i)igN being an independent family of real valued Brownian motions. The correlation function 
of the process W is then given by 

E{W{t,x)W{s,y)) = c{x -y){s At), x,yeR, s,t>0, 

where 

c{z) = / k{z + u)k{u)du. 
Jr 

The existence and uniqueness of solutions for stochastic KdV equations of the type (jl.ip 
but with an additive noise have been studied in [3], [7], [8]. The multiplicative case with 
homogeneous noise as described above was considered in [6]: assuming, together with the 
above condition, that k is an integrable function of a; G M allowed us to prove the global 
existence and uniqueness of solutions to equation (jl.ip in the energy space (M) , that is in 
the space where both the mass 

(1.4) m(n) = - I u^{x)dx 

2 Jr 

and the energy 

(1.5) H{u) = 11 {d^ufdx -If u'^dx 

are well defined. Note that ra and H are conserved for the equation without noise, that is 

(1.6) dtu + dlu + ]^d^{u'^) = Q. 
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Under the above conditions on k, it was then proved in [6j that for any given initial data 
uq G H^CM.), there is a unique solution u of (jl.ip with paths a.s. continuous for t € with 
values in H^{R). 

Our aim in this article is to analyze the qualitative influence of a noise on a soliton solution 
of the deterministic equation. More precisely, we study the qualitative behavior of solutions 
of (jl.ip in the limit e tends to zero, assuming that the initial state of the solution is a soliton 
of equation (|1.6p . We recall indeed that equation ()1.6p possesses a two-parameter family of 
solitary waves (or soliton) solutions, propagating with a constant velocity c > 0, with the 
expression iic,xo(^) ^) = ^c{x — ct + xq), xq G M, where 

(1.7) ^.(x) 



2cosh^(VcT 



satisfies the equation 



(1.8) ^'^ - + \pI = 0. 

We do not recall here the well-known results concerning the stability of the soliton solutions 
Uc,a;o ill equation (jl.6p , but we refer to [2] , [3] , or ^18j for a review of the stability questions 
using PDE methods, or to [13] and [20] for a review of the stability of the solitons with the 
help of the inverse scattering transform. 

Let us consider as in [5j the solution u'^(t,x) of equation (jl.ip which is such that ii^(0, x) = 
(po{x) where cq > is fixed. Then, in Section 2, we show, as we did in [5j for the additive 
equation that up to times Ce~^, where (7 is a constant, we may write the solution as 

(1.9) u%t,x) = ip,.(^t^{x-x%t)) + er]%x-x%t)) 

where the modulation parameters d^{t) and x^{t) satisfy a system of stochastic differential 
equations and the remaining term er]^ is small in H^{M.). We then prove in Section 3 that the 
process ry^ converges as e goes to zero, in quadratic mean, to a centered Gaussian process rj 
which satisfies an additively driven linear equation, with a conservative deterministic part; we 
also investigate the behavior of the process ry as t goes to infinity and prove that rj is in some 
sense a Ornstein-Uhlenbeck process, with a unique Gaussian invariant measure. In addition, 
the parameters x^{t) and c^{t) may be developed up to order one in e and we get 



J dx^ = codt + eBidt + edB[ 
\ dd" = edBi +o{e), 



■2 + o{e) 



where Bi and B2 are correlated real valued Brownian motions; keeping only the order one 
terms in those modulation parameters, we then obtain a diffusion result on the modulated 
soliton similar to the result obtained by Wadati in [21j, but with a different time exponent (see 
Section 4). 

In all what follows, (., .) will denote the inner product in L^( 



{u,v) = / u{x)v{x)dx 
Jr 

and we denote by T^q the translation operator defined for if £ C(M) by {Txgif){x) = if{x + xo). 
Note that since the process W is stationary in space, for any xq G M the process T^qW is still 



4 



A. DE BOUARD AND A. DEBUSSCHE 



a Wiener process with covariance <j)(j)*. Indeed by ()1.3p . 



%,W{t,x) = Y,iMix + xo)(3k{t) = Y,i^^k){x)(3kit) 




with ek{x) = TigCfc. 



2. Modulation and estimate on the exit time 



In this section, we prove the following theorem. 

Theorem 2.1. Assume that the kernel k of the noise satisfies together with k £ L^{]R.) 

and let cq be fixed. For e > 0, let u^{t,x), as defined above, be the solution of with 
it(0,a;) = lPcq{x). Then there exists oq > such that, for each a, < a < a^, there is a 
stopping time > a.s. and there are semi-martingale processes c^{t) and x^{t), defined a.s. 
fort < T^, with values respectively inM"*"* andM, so that if we set£r]^{t) = u'^{t, .+a;^(f))— 
then a.s. for t < t^, ||e?7^(t)||i < a and \c^{t) — co\ < a. In addition, for oq sufficiently small, 
and any a < uq, there is a constant C > 0, depending only on a and cq, such that for any 
T > 0, there is an Eq > 0, with, for each e < eq, 



It was noticed heuristically in [5j, and proved in [lOj that in the additive case, the use of 
the modulation parameters x^{t) and c^{t) was necessary in order to get the estimate (|2.ip . 
Indeed, it was proved in [TO] that if we denote by fa" = inf{t > 0, .) — </?colli > ck}) 

where u^'" is here the solution of equation (|l.ip . but with an additive noise that becomes 
stationary in space as n goes to infinity (see [10] for a precise statement) then there exists a 
constant C(a,co) which depends on a and cq but not on T such that 



It is not clear that (j2.2p is still true in the present multiplicative case, because the proof involves 
a controlability problem with a potential which - up to now - is open. 

Note also that the decomposition given in Theorem 12. II is not unique, and is determined by 
the choice of specific orthogonality conditions (see the proof below). In particular, contrary to 
the additive case, we will be able here to investigate the asymptotic behavior in time of the 
limit process by choosing one particular decomposition of the form given in Theorem 12.11 This 
is the object of Section 3.3. 

Proof of Theorem \2. il The proof follows closely the proof of Theorem 2.1 in [5j and we refer to [5] 
for more details. The parameters x^(t) and c^{t) are obtained thanks to the use of the implicit 
function Theorem. These are then local semi- martingales defined as long as |c^(t) — cq\ < a 
and ||M^(t, . + x'^it)) — ^coWi < and setting 



(2.1) 




(2.2) 



lim„_.^lim,^o^^' logP (f^'^ < T) > 



C(a,co) 

r3 



e7]'{t) = u'{t,. + x'{t))-ip,e^t) 




(?7^</'co) = (^/^5x'/3co) = 0. 



SOLITON DYNAMICS FOR MULTIPLICATIVE KDV EQUATION 5 

In order to estimate the exit time 

= inf{t > 0, \d^{t) - Col > a or ||er?'^(t)||i > a}, 
we make use , as in [5], of the functional defined for u G H^{M.), 

(2.4) Qcoiu) := H{u) +com{u) 

where H and m are defined respectively in (jl.4p and (jl.Sp . Note that (pcg is a critical point of 
QcQ- We denote by Lc^ the linearized operator around ipc^, that is 

(2.5) L,, = -dl + co-2ip,,. 

The next lemma, which is proved with the use of the Ito Formula, using the same regularization 
procedure as in [4J, gives the evolution of H and m for the solution of (jl.ip with m^(0) = ipcg ■ 

Lemma 2.2. For any stopping time t < +cxd a.s, one has 

rr rr 

m{u'iT)) =m{ipco) -e {iu'-){s),dW{s)) + e^\k\l2 m{u%s))ds 

Jo Jo 

and 

(2.6) Hiu%T)) = H{^^) + e r {d^u',d^[u'dW[s)))-'- [\{u')\dW{s)) 

Jo ^ Jo 

(2.7) +^£||fc|2^|a,n-|2^ + |A;'|2,|^^|2^|ds 

(2.8) -^ly^ r f {u'f\cpek\'dxds. 

^ Jo Jr 



Consider v > such that {Q'^^{ipco)v,v) > lyWvWf for any v G satisfying (f,(^co) = 
{v, dx^co) = 0- The existence of such a constant is a classical result (see [2] or [3j). Then it is 
easy to show (see [5]) that there is a constant C(ao) > such that for any t < t^, 

(2.9) Qcoiu'it, . + x^m - Qcoi^,.^t)) > jlWml - C\c%t) - cop. 
Now, if r = A t, then by (j2.9p . the translation invariance of Qcq, and Lemma 12.21 

\\ev'{r)\\l < - [QcoM-Qcoi^c^ir))] +e [ {d^u' {s) , dx{u' dW {s))) 

^ T 2 r 

r{{u^f{s),dw{s))+'- n\k\udxu^\i,+\k'\uu^\i,)ds 

(2.10) y° ... ^ ^° 

-tE/ {unHsMek\^dxds-coe {{u'f,dWis)) 
^ f, Jo Jr Jo 

+coe''\k\l, r m{u^s))ds + C|c"(r) - cop. 

The term |c^(t) — co| is then estimated thanks to the orthogonality condition (r]^,(pcg) = 
and the evolution of m{u^{T)) given in Lemma [2^ one obtains, for some constants ^ > and 
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C > 0, depending only on cq and oq (with a < oq) 



file" (t) - Co\ < 



\v>co\h 



\fcHT)\L2\ 



< 



\^r]%T)\l2 + Calc^ir) - Col + 2e 
+2e^\k\l. 



iiu'y,dWis)) 



\u^(s)\l2ds. 



Hence, choosing ao sufficiently small one gets 



(2.11) 



\c%t)-co\' < 
which, once inserted into (j2.10p leads to 



\eri%r) 



< C 



\u^is)\i2ds 



id,u',d,iu'dWis))) 



{{u'r,dWis)) 



+ CqE 



{{u^)^dWis)) 



+e^\k\ 



+ e^\\kfi 



L2 



{(un^dwis)) 

\\u'{s)\\lds + e^\k\l. 



lids 



\u''{s)\l2ds 



With this estimate in hand, together with (j2.1ip . the conclusion of Theorem 12.11 follows with 
the same arguments as in the proof of Proposition 3.1 in [lOj. These arguments rely on classical 
exponential tail estimates for stochastic integrals, after noticing that ||?/^(s)||i < C, a.s. for 
s S [0, A T] and a < oq, so that the quadratic variation of each of the integrals involved in 
the above estimates are bounded above by CT. □ 

3. A CENTRAL LIMIT THEOREM 

This section is devoted to the proof of the next theorem: 

Theorem 3.1. Under the assumptions of Theorem \2.1\ let a < ao be fixed. Then we can 
find (?{t) and x^{t) satisfying the conclusion of Theorem \2.1\ such that if fj^ is defined as in 
Theorem \2. 11 for any T > 0, the process {fj^ {t))te[o,T] converges in L?'{Q.] L°°(0, AT; L^(M))) 
to a Gaussian process fj satisfying the additive linear equation 

(3.1) df] = dxLcafjdt + QifcodW, 

with fj{0) = 0, where W is the Wiener process with covariance (jxp* given by W = Tc^tW , and 
Q is a projection operator. Moreover, for a > sufficiently small compared to cq, the process 
w{t,x) = e°'^fi{t,x) is a well defined valued process, of Ornstein- Uhlenbeck type, which 
converges in law to an -valued Gaussian random variable as t goes to infinity. 

The conclusion of Theorem 13.11 will be obtained in three steps. The first step consists in 
estimating the modulation parameters obtained in Theorem 12.11 in terms of ry"^, using the 
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equations for those parameters; then the convergence of r]^ as e tends to zero is proved, and 
finally in the third step, a slight change in the modulation parameters is performed, in order 
that the limit process 77 may be written as an Ornstein-Uhlenbeck process. 

From now on, we assume that a is fixed and sufficiently small, so that the conclusion of 
Theorem 12.11 holds . and we denote by r^. 

3.1. Modulation equations. Since we know that the modulation parameters x^{t) and c^{t) 
are semi-martingale processes adapted to the filtration generated by {W{t))t>o, we may a priori 
write the stochastic evolution equations for those parameters in the form 

, . j dx^ = c^dt + ey^dt + dW) 

^ ' \ dc^ = ea^dt + e{b^,dW) 

where and are real valued adapted processes with a.s. locally integrable paths on [0, r^), 
and 6^, are predictable processes with paths a.s. in Lf^^{0,T^; L'^ We then proceed as 
in [5j : the Ito-Wentzell Formula applied to u^{t,x + together with equation (jl.ip for 

and the first equation of (j3.2p for x^ give a stochastic evolution equation for u^{t,x -\- x"^). 
On the other hand, the standard Ito Formula together with the second equation of (j3.2p for 
give an equation for the evolution of (/?c=(j). Replacing then V9cs(t) +£r]^(t,x) for u'^{t,x + x^{t)) 
in the first equation leads to the following stochastic equation for the evolution of ?7^(t) : 



d-i]"^ = dxLcoV'^dt + (y^dx^Pc^ - a'^dc^c^)dt - dx{{^c^ - ipco)'n^)dt 
+(c^ - CO + ey')dxr]'dt - |9,((r/^)2)dt + ipc^%edW 
^3 3^ +dx^c^iz', dW) - dcipc^b', dW) + £7]'%.dW + ea,r?^(z^ dW) 

+^dlipcA(t>* z'\l2dt - f 9,VcH'A*^1i2dt + e ^ dx{^c-Txe(t>ei){z' , (^ei)dt 

^e^dlr^'\(t>*z%,dt + e^Y. dM%^^ei){z', ^ei)dt 

where Lcq is defined in (|2.5p . Now, taking the L^- inner product of equation (j3.3p with ipc^, on 
the one hand, and with dx'-Pco on the other hand, then using the orthogonality conditions (j2.3p 
and the fact that LcqOx^Pcq = 0; finally identifying the drift parts and the martingale parts 
of each of the resulting equations lead to the same kind of system that we previously obtained 
in [5]; namely, setting 

and Zm-(f^:^ 

then one gets for the drift parts 

(3.4) A%t)Y%t) = G%t) 

where 



(3.5) A%t) 



(dx^Pc^ + edxTj", dx^Pco) -{dc'Pc^,dx^cQ, 

- {dx^c'^ , (fco ) {dc<fc^ , <fco ) 
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and 



with 



Glit) = (7?^Leoa2vPeo) + (c^-co)(r?^^2<^co) + i(9.(r?^)^9.(^c„) 



and 



G\{t) = -|(5x(??^)^¥'co) - {dxii^e^ - 'Pco)ri^),Vco) + i{^x<Pc^,<PcoM*Z^\l2 

(3.7) -U9'^'Pc'^,^co)\4>*b''\l2 + £^{z',(t)ei){dx{^Pc^%'^4>ei),ipco) 

note that A%t) = + 0(|c^ - co| + He??"!!!), a.s. for t < with 

V i^co,dc'Pco)J 

and 0(1 — Col + 11??^ Ill) is uniform in and to as long as t <t^. Concerning the martingale 
parts, one gets the equation 

(3.8) A%t)Zf{t) = Ff{t), V/gN 
with 

(3.9) F^{t) = H^"^^' + ^rj2%^fl, dx^co)\ 
^ ^ V Wc^ +£rj^)%^(f>ei,<Pco)- J 

Proposition 3.2. Under the above assumptions, there is a constant ai > 0, such that if 
a < ai, then 

(3.10) \(l)*z%t)\L2 + 10*6^1^2 < Ci\k\L2, a.s. for t < 
and 

(3.11) \a%t)\ + \y%t)\ < C2\ri'{t)\L2 + eCg, a.s. for t < t' 

for some constants C\, C2, C3, depending only on a and cq, and for any £ < Eq. 
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Proof The proof is exactly the same as the proof of Corohary 4.3 in [5j, once noticed that, a.s. 
for t < r^ 



< C V / {if^e + e7i'f{x)[{%.k) * ei]\x)dx 

I -^^ 

< / {<fc^ + eri'fix) y2{%.k{x - .), eifdx 
Jr I 

< C [ {^c^+eri')\x)\%.k{x - .)\l,dx 

< C\k\U^,.+erj%,<C\k\l, 

where we have used the Parseval equaUty in the fourth Hne. □ 

3.2. Convergence of r/*". Let us first assume that r/*" has a hmit as e goes to zero, and take 
formahy the hmit as e goes to zero in the preceding equations. Then, as was noticed above. 



hence 



1 



(3-12) hni0*2;^ = -— -^{Tc,,t(l)y {fcodxVco) ■= z 

(3.13) hmr^^ = 7 1 -{%,tcl^*)ivl) := b 



(3-14) Hm/ = — -^{'n,Lcodlipco) := y 

and 



1 

Il2 



(3.15) hm = 0. 
Moreover, formally, r] satisfies the equation 

drj = a^Lco??dt + |^-^i-|7^(7?,Lco9^v3co)9x'/'codi 

(3.16) +ipcoTcotdW - 2\dJ,^\l^ {dx{^l^),%otdW)d,x^co 

- (^.oi^.,M o ' '^cotdW)d,ip,, . 

It is easy to show that (j3.16p has a unique adapted solution t] with paths a.s. in C{M.'^ , H^) 
satisfying 77(0) = 0. Moreover using the fact that {dc^coi dx^co) = 0, one easily gets from the 
above equation that (rj, ip^o) = (rj, dx^co) = 0, Vi > 0. 

Next, we make use of the following lemmas, whose proofs are obtained in the same way as 
the corresponding Lemmas in [5]. 
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Lemma 3.3. Let r] be the solution of 1^3. 16\) with r/(0) = 0. Then, for any T > 0, there is a 
constant C depending only on cq, T and \\k\\i such that 

E{ut)\\i)<c, yt<T. 

Lemma 3.4. Let ry^ be the solution of Ii3.3\) . defined for t G [0,t'^[, obtained thanks to the 
modulation procedure of Section 2. Then, for any T > 0, 

e( sup W(.t)\l2) <C(r,Q,co,||A:||i). 

^t<T'AT ' 

The above lemmas show that 

(3.17) VT > 0, Vg > 2, hmEf sup |c^(t)-co|'^ 

Indeed, the expression of c^{t) — cq given by ()3.2p together with (IS.lOp and (jS.lip imply easily 

e( sup \c%t) - co\^) <Ce^[l+E \r]%s)\l2ds] 

with C = C(a, Co, T, || Then, (j3.17p is deduced form Lemma [3^ for q = 2, and follows for 
all other values of q from the uniform boundedness of |c^(t) — cq\ on [0, T A r"^]. Note that an 
immediate consequence of (|3.17p is the fact that 

(3.18) VT > 0, Vg > 2, limE^ sup W'PcHt) - "PcoWt] = 0- 

We will finally need the next lemma. 
Lemma 3.5. For any T > 0, and any q>l, 

limE( sup (Y,\Zt{t)-Zi{t)\' 

where we have set for I G N 



■t<TAr^ ZgN 



z and h being given by i3.12\) and i3.13\) . respectively. 

Proof Here again, it is sufficient to consider the case q = 1- We recall that Zf satisfies equation 
(j3.8p . First, it is clear that 



limEf sup \\{A^t))-^ - (Ao(t))"i|p9) = 0, > 1. 

On the other hand, in view of (j3.9p . denoting F^(t) the formal limit of F^i^t), one has 
e( sup 5;|Ff(t)-F0(t)p) 

<CE( sup ^\dx(Pco{%^(t> -%ot(t>)ei\'i2j + CE( sup \\iPc^{t) - fcoWi 



t<TAT^ ^ ' ^ t<TAi 
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and 

t<TAT^ 



E( sup S2\dx(pcoi%^(t> -%ot<P)ei\l2 
< ||(PcoIIiIe( sup \k{. + x\t)-CQt)-kf^^. 



■ t<TAT^ 

Then, the Ito Formula apphed to the function 

/C^(i, x) = {k{x + x'it) - Cot) - k{x)f 

using equation (j3.2p for dx^{t), together with (jS.lOp . (|3.1ip . and (j3.17p lead to the conclusion 
of Lemma 13.51 □ 
Now, in order to prove that 

(3.19) limE( sup \7]'{t)-r](t)\l,) =0, 

where rj is the solution of (|3.16p with 77(0) = 0, it suffices to set = ij^ — ij, to deduce from 
(j3.16p and (j3.3p the equation for dv^ and to apply the Ito Formula to get the evolution of 
|^^^|^2- We do not give the details of those tedious, but easy computations. Finally, the use of 
the following estimates : 



e|(^;^a,((^?^)2))| =e|(a,.r?,(7?^)2)| <e||7?||i|r7^|2 



< Ce\\r|\\l\rJ'\%^\^^l^'\%^ < C^\\ri\\i\r,'\%^ 



on the one hand, and 

-y| + |a^| < C{\v''\l2 + - coM\l2 +e|r?^|^2 + W\lA\Vc'^ - "fcoh+e) 

which is obtained as in the proof of Lemma [3. 5 1 on the other hand, together with Lemma [3. 3 1 to 
l3.5l allow to get the conclusion, that is the convergence of r/*" to r/ in Lp'iQ., L°°(0, AT; ^^(M))). 
□ 

3.3. Complements on the limit equation. First of all, we note that the modulation equa- 
tions may be written at order one in e as 

dx" = codt + eydt + eWidt + edW2 + o(e) 
dc" = edWi + o{e) 

where 

y = \d^(pco\l2iv,Lcod'^(Pco), 
Wi{t) = {ip,„d,^,,)-\vl,W{t)) 

and 



1 

'2' 

Note that Wi and W2 are real valued Brownian motions, which are independent since 



nwi{t)W2is)) = -^\d,^co\lH^co,dc^co)~\4>*idAvl)),<P*i^^^^^^^^ 

because the operator (p* commutes with spatial derivation. 

Now, we want to investigate the asymptotic behavior in time of the process rj. However, in 
the present form, the process rj does not converge in law as t goes to infinity; this is due to 
the fact that the preceding modulation does not exactly correspond to the projection of the 
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solution on the (two-dimensional) center manifold, in which case the remaining term would 
belong to the stable manifold around the soliton trajectory. We now show that by slightly 
changing the modulation parameters, we can get a new decomposition of the solution which 
is defined on the same time interval as before, but which fits with the preceding requirements. 
For that purpose, we first need to recall a few facts from |18j . 

The generalized nullspace of the operator dxLcQ (that is the operator arising in the linearized 
evolution equation in the soliton reference frame) is spanned by the functions dx^co ^^id ^c'-Pco, 
with the equality 

dxLcodc^Pco = -dxVco 
and there are constants 6i and 62 (with 61 = {ipcQ,dc^Pco)) such that if we set 

/X 
dc(pco{y)dy + 92(pco and g2{x) = OupcQ 
'OO 

then the generalized nullspace of —L^^dx is spanned by gi and ^2 and 

{gi,dx^co) = 1, (5i,9c(/7co) = 0, {g2,dxVco) = 0, {g2,dcfco) = 1- 

We also set, for a > 0, 

(x) = e'^^dx^co, f2{x) = e-^d.ip,,, g1{x) = e~'^^~gi{x), gl{x) = e-^^U^), 

so that (/", g'j) = 5ij. Then the operator Aa defined for a > by Aa = e°-^dxLcf,e~"'^ has a well 
defined generalized nullspace spanned by ff, /I" and the spectral projection on this nullspace 
is given by Pw = X]fc=i(^j 5^)/^ where w = e'^^v, and v is an function. Moreover, if 
Q = I — P, then Q is the spectral projection on the stable manifold of Aa, and under the 
condition < a < y^co/S, there are constants C > and 6 > such that 

(3.20) ||e^"*Qu;||i < Ce~^*||u;||i, Vt > 0, Vu; E H^, 

where e"^"* is the C'^-semi-group generated by Aa (see Theorem 4.2 in |18j). 

Now, let r] be the solution of ()3.16p with r]{0) = 0, and consider w{t,x) = e"'^rj{t,x). Note 
that the orthogonality condition (77, (/9co) = implies (1^,52) = 0, so that Pw = X{t)ff with 
X{t) = {w{t),gi) a real valued stochastic process whose evolution is given by 



(3.21) 



Kt) = I \dxtPco\J{v{s),Lc^,dlipco)ds- \dxiPco\Ji^codxiPco,dW{s)) 



{e''^ip,,dW{s),g1) 



where we have used (|3.16p and the fact that AaPw = and A(0) = 0. Hence, X{t) is bounded 
in L'^(17;L°°(0,r Ar^)) by LemmaEH Let us set x%t) = x%t)-eX{t) for t E [0,r^[. Then 

(3.22) u%t,x + x%t)) = ipc-^ij^^ix) + e7)''{t,x) 

with 

fi^{t,x) = -{ifc'^(^t){x -e\{t)) - ipce{t){x)) +r]%t,x - eX{t)). 
Note that, a.s. for t < : 

\(Pc^(t){- - eA(t)) - ipce(t) - £:Kt)dx^c^(t)\L^ < e'^X^{t)C{co,a). 
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Hence, it follows from Lemma 13.31 13-41 and the above bound on A that 

(3.23) limE( sup \f{t)-fi{t)\l2)=0 

with f/(t) = r/(t) — X(t)dx^co- So now, with this new decomposition, we clearly have, setting 
w{t, x) = e'^^fj^t, x) : 

Pw = 0, Qw = Qw. 
Also, if W2 = Qw, then the equation (j3.16p implies 

(3.24) dw2 = AaW2dt + Qe'^^'ipcodW 
hence 

Jo 

the trace of the covariance operator of the Gaussian process W2 in may be easily computed 
and estimated thanks to ()3.20p as 



r V lle^-'^Qe'^^^co'/'e/llfdCT < C( f e'^" da) V 
Jo , Jo , 



_ e'^'^Vco^eillida < C\\k\\i\\e''^^co\\i- 

I 

Moreover, this covariance operator converges as t goes to infinity and it follows that W2 con- 
verges in law in to a Gaussian random variable. The end of the statement of Theorem 13.11 
follows, setting Qv = e~°-^Qe°-^v. □ 

4. A REMARK ON THE SOLITON DIFFUSION 

Let us go back to the stochastic evolution equations for the new modulation parameters, 
that we may write as 



(4.1) 



dx^ = codt + eBidt + edB2 + o(e) 
dc^ = edBi + o{e) 



with Bi = Wi and B2 = -(e"^93co^^(0. 5i ) = -(^^(*), V^coffi)- Note that Bi and B2 are now 
correlated Brownian motions. We denote by 

cr = {(7ij)ij = cov{Bi,B2). 

If we keep only the order one terms in e i.e. we consider the solution (X^(t),C'^{t)) of the 
system of SDEs 

dX^ =codt + eBidt + edB2 



dC^ = edBi, 

then {X^{t) — cot,C^{t) — cq) is a centered Gaussian vector, and it is easy to compute its 
covariance matrix. Let us denote by /if the law of (X^(t) — cot, C^{t) — cq); we may compute 



maxE((/^c-(t)(x-X^(t)) 



(4.2) / / ^c+co(2; - Co* - y)/if(dy,dc) 



max . . / / ^c+co{x - Cot - y) exp - ^ f ^ j . (f \ ^dcdy 



xei (dets)V2 7 7 -'---V 2" \yj-\y. 
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where S is the covariance matrix of {X^{t) — cot, C^{t) — cq), given by 

/ \ 
IP ^2 I O'll* o"i2i + o'ii-T 1 

^ = e I t-2 ^ .3 



t 2 t 

Cl2i + Clly (722t + 0"i2t +'7ll3". 

It is not difficult to see that 

Inserting this inequality in ()4.2p . using the fact that '^c{x) = c</9i(a/cx) and integrating in y 
give the bound 

where is a constant, and since 



I-+00 ^2 

/ ^/ce ^dc< Ka^'^ 
Jo 



for another constant K, it follows 



(4.3) maxE((^c;.(t)(x - < Koe-'^/^t-^^^ 



for t large enough. 

This inequality has to be compared to the result of [21j where an additive equation with a 
white noise in time was considered. An inequality of the form (j4.3p was obtained, but with a 
power t"'^/^ instead of t~^^'^. 
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